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ABSTRACT 

We present a large collection of photometric data on the Blazar PKS 0537-441 in the VRIJHK 
bands taken in 2004—2009. At least three flare-like episodes with months duration, and >3 mag 
amplitude are apparent. The spectral energy distribution is consistent with a power law, and no 
indication of a thermal component is found. We searched for short time scale variability, and an 
interesting event was identified in the J band, with a duration of ^25 minutes. 

Subject headings: galaxies: active — galaxies: BL Lacertae objects: PKS0537-441 



1. Introduction 

Blazars are active galactic nuclei (AGN) dom- 
inated by a relativistic jet pointing to the ob- 
server. The relativistic beaming enhances the con- 
tinuum, and therefore the equivalent widths of 
spectral lines produced by photoionization of cold 
gas clouds around the nucleus may be modest, 
comparing with other AGN classes. In BL Lac ob- 
jects, line weakness is extreme (equivalent width 
<5A). 

In order to relate dimensions and time scales 
in the jet comoving frame to what is directly ob- 
served, one should go through Lorentz transfor- 
mations, which are characterized by the Doppler 
factor defined as 5 = r _1 (l — /3cos0) -1 , T is the 
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Lorentz factor and 9 is the angle of sight. When 
the relevant dimensions are much smaller than the 
spatial angular resolution, which is the usual case 
for the inner jet, the properties of the jet can be 
effectively constrained by the study of the va ri- 



ability of the emission (e.g. lUlrich et al. II1997I ). 



Blazars are generally observed in a large fre- 
quency range, which, in some cases, covers the 
entire electromagnetic spectrum from the radio 
band to the TeV gamma rays (e.g. Fossati et 
al. 1998, Donato et al. 2001, Ghisellini et al. 
2010). The basic model for explaining the spec- 
trum of blazars is in terms of synchrotron radi- 
ation, plus an inverse- Compton (IC) component, 
which may be due to scattering off the electrons 
of the synchrotron photons, or those coming from 
the accretion disk, broad-line region or infrared 
torus. In practice, spectra are well constrained if 
the source size can be obtained though variability 
studies. The ideal situation would be to obtain 
the time dependence of the spectral energy distri- 
bution in the entire band of observation, and this 
has been done for very few objects (e.g. Hartmann 
et al. 2001, Ballo et al. 2002, Acciari et al. 2009, 
Fuhrmann et al. 2010). 

The systematic study of blazar variability in the 
optical and infrared bands has greatly profited of 
the development of robotic telescopes, which en- 
sure long observing campaigns, with relatively lit- 
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tie intervention from the observer. The large ob- 
serving time can yield information complement- 
ing that of multifrequency campaigns, which ob- 
viously can cover at most few days. 

Here we report on 2004—2009 six filters pho- 
tometry of the blazar PKS 0537-441, which is un- 
precedented for the observing time dedicated to 
the source. PKS 0537-441 is a z = 0.896 bright 
(V~13) blazar, which was studied from radio to 
GeV 7 rays. The main characteristics of the 
broad-band spectrum are the synchrotron peak 
in the IR, the IC hump with a maximum in the 
7 rays, the X-ray spectrum sampling the tail of 
the synchrotron and the rise of the IC, and the 
presence of emission lines in the optical spectrum. 
Therefore, the object belongs to the low-frequency 
peaked subclass of blazars (Pian et al. 2002, 2007). 

In Section 2, we describe the observational pro- 
cedures, data reduction and analysis. The results 
are given distinguishing between long- and short- 
term variability (Sections 3 and 4). In the discus- 
sion section in particular we focus on the relevance 
of short time scales of variability in constraining 
the beaming of the source. 

2. Telescope and data analysis 

The Rapid Eye Mount (REM) telescope is a 
60 cm fully robotic instrument located at the ESO 
Cerro La Silla observatory (Chile). The telescope 
has a Ritchey-Chretien configuration with a 60 
cm f /2.2 primary and an overall f /8 focal ra- 
tio in a fast moving alt- azimuth mount providing 
two stable Nasmyth focal stations. At one of the 
two foci, the telescope simultaneously feeds, by 
means of a dichroic, two ca meras: REMIR f or the 
Near-InfraRed (NIR) (see IConconi et al. 1 12004 



and ROSS (see iTosti et al. 



120041 ) for the optical. 
Both the cameras have a field of view of 10' x 10' 
and imaging capabilities with the usual NIR (z 1 , J, 
H, and K) and Johnson-Cousins VRI filters. The 
telescope was originally designed for follow-up of 
Swift gamma-ray bursts, but for most of the time 
it is used for standard observational programs. 
More information about the REM project and ca- 
pabilities can be found in Zerbi et al. (2001), Chin- 
carini et al. (2003), Covino et al. (2004). 

Both REMIR and ROSS instruments were used 
in order to obtain nearly simultaneous data and 
to study the spectral behavior of PKS0537-441 at 



different levels of flux. The REM telescope has 
systematically observed this source for 5 years be- 
tween December 2004 and March 2009. 

Data analysis has been carried out following 
standard recipes. Instrumental magnitudes were 
obtained via aperture photometry by means of 
an automatic pipeline^ which aligns and rotates 
frames to a common reference, derives and cali- 
brates photometry comparing to standard stars in 
the field from the 2MASS catal oguH in the NIR 



and fr om sequences reported by lHamuv &: Maza 



(|1989i ) in the optical. Results were controlled man- 
ually using the Graphical Astronomy and Image 
Analysis (GAIAf] tool: we selected only frames 
with reliable photometry and absolute calibration 
based on the presence of the calibration stars, 
on their consistency, stable background, correct 
centering, etc. In particular, for all frames, we 
required that at least three out four calibration 
stars were detected. The reference stars, the same 
in the optical and in the NIR, are stars Nr. 1, 

2, 3 and 6 as labeled in Hamuy & Maza (1989) 
corresponding to the 2MASS stars 05390451 - 
4406382, 05384933-4401300, 05385997-4401191 
and 05383780 - 4403241, respectively. The agree- 
ment between the automatic pipeline and the 
manual analysis is typically better than 10% in 
flux. 

In the following, fluxes are corrected with Ay = 
0.18 (Kalberla et al. 2005) and the extinction laws 
of Cardelli et al. (1989), while magnitudes are the 
observed ones. The quoted errors are at la and 
they include photometric and calibration uncer- 
tainties. 

3. Long time scale variability 

Fig. Q] represents the light curves of the source 
in the six bands. Part of these data was already 
published in Dolcini et al. (2005), Pian et al. 
(2007), Impiombato et al. (2008). Pucella et al. 
(2010) consider some photometry of the source in 
October 2008 coincident with an AGILE 7— ray 
observation. The 2006 data, which are noticeable 
for the large and complex variability, are reported 
in an enlarged version in Fig. [2] For all photo- 



1 http: / / www.merate.mi.astro.it / utenti / covino / usermanual.html 
2 http:/ /irsa. ipac.caltech.edu 

^http: / /star- www. dur.ac.uk/ ~pdraper/gaia/gaia.html 
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1— Optical and NIR light curves of PKS0537-441 in the 2004-2009. 
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metric points (see Table [T]), we give the filter, the 
reference MJD at the beginning of the exposure, 
the integration time, the magnitude, and the cor- 
responding error. The typical integration times 
are 20—600 s in the optical and 170 s in infrared 
bands. 

The maximum and minimum states are re- 
ported for each filter in Tabled The values have 
been controlled through the manual analysis and 
are in good agreement with the automatic analysis 
ones. The variability ranges from ~ 4 mag in V 
to ~ 6 mag in J. Note that the time coverages in 
the various filters are not equal. 

4. Short time scale variability 

Since the typical exposure duration is 20—600 
seconds, variability on a few minute time scale 
can be searched for. However, our observing 
campaigns were not optimized for the study of 
such variability and, in fact, there are relatively 
few epochs where 5—10 points within an hour 
are available. Moreover, the observing procedure 
favours filter shifts after few photometric points. 
Nonetheless, we have considered the usual defini- 
tion of an exponential increase/decrease of flux: 



F(t) = F(t )exp[-(t-t )M 



(1) 



where F(t) and F(to) are the fluxes at the time 
t and to, respectively, and r is the characteristic 
time scale. The time scale r intrinsic to the source 
is then calculated as r- mt — r/(l + z). 

We considered only frames with reliable pho- 
tometry basing on consistency among the various 
reference stars (at least three) and no image de- 
fects close to the target position. An automatic 
procedure was launched in order to find events of 
minimal time scale, r (< 1 hour), with flux differ- 
ences with a significance above 5cr. One event was 
recovered, which appear of interest. It is reported 
in Table [31 all the frames were carefully checked 
manually. No simultaneous events satisfying our 
criteria were detected in the optical bands. Note 
that the telescope is designed to provide the best 
sensitivity in the NIR. 

5. Discussion 

The light curves of 2004-2009 (Figs. ED and 
[2]) correspond to rather irregular exposures and 



are interrupted by large gaps, partly due to the 
visibility conditions of the source. Nonetheless, 
this is the largest collection of published data on 
PKS 0537-441 in terms of overall exposure (t ~ 
3 x 10 5 s) and number of photometric points in all 
bands (about 4000 points). 

A visual inspection of Figs. [T] and [2] indicates 
very structured light curves, with grossly similar 
shapes in all bands. A number of rapid bright- 
enings and fadings, which may suggest flare like 
episodes, are apparent. Typical durations are 
months and the amplitudes are of several mag- 
nit udes (see also Table [2t ■ The event de s cribed 



by iDolcini et al~l (|2005l ) and IPian et al. I (|2007f ) 



is that of Jan 2005 (days 41-71 in VRI bands, 
see Fig. [[]). That of February-March 2008 (days 
1189-1216 in VRI JHK bands, see Fig ED) is the 
one considered by llmpiombato et al (|2008l ) . A 
double flare structure is present in the period 
21/02/2006-31/12/2006 (see F ig. H) Proba bly 
the event of 1971, reported by Eggenl (Il973h . 4 



mag in V in 8 months, is analogous to those ex- 
amined here. 

We compared the spectral shapes at three 
epochs where all six filters are available and the 
intensity states are significantly different. An in- 
teresting case is shown in Fig. [3l where it is 
apparent that the SED can be satisfactorily de- 
scribed by power laws. It is noticeable that the 
higher and lower states, which are separated by 
two weeks, exhibit an intensity variation by a fac- 
tor ~ 10 and the spectral index changed only 
slightly (Aa ~ 0.1). Therefore, in both the low 
and high state the emission appears dominated by 
the jet, and no evidence of a thermal component is 
apparent. Some deviations from a power law may 
appe ar in the low-state 2004 December 21 (Fig. 



3 of Dolcini et al. 



2005h. The ther mal compo- 



Pian et al. I (|2007h is probably 



nent mentioned by 
a consequence of the extreme variability of the 
source. 

We now concentrate on the short time scale 
variability, and in particular on the event of 2006 
October 30 (MJD 54038). This is rather rapid: 
comparable to what observed by Villata et al. 
(2008) in the notorious highly variable BL Lac ob- 
ject S5 0716+714, about 0.002 mag/min in V. A 
time scale T ln t ~ 25 min corresponds to a size 
in the source frame R' < CT lnt ~ 4 x 10 14 cm, 
having considered S ~ 10, typical of blazars (e.g. 
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Fig. 2.— Optical and NIR light curves of PKS0537-441 in the year 2006. 



■5 




6 



Ghisellini et al 2010, and references therein). This 
can be compared with other constraints on the 
si ze of the source. I ndeed, as shown in Table 4 
of iPian et al. I ([20051 ) . the mass estimates of PKS 
0537-441 range from 0.5 x 10 8 M Q to 16 x 1O 8 M . 
They are based on the measurements of line widths 
and the continuum intensity, which enable one 
to deduce the velocity of the clouds emitting the 
broad lines and their distance. The mass then fol- 
lows from the virial theorem. The problem is that 
the continuum is dominated by the jet and one 
should use relativistically corrected values. Such 
a correction is introd uced, but in a somew hat gross 



formulation, only bv IPian et al. (20051), who fa- 



vors the lowest value of mass. The gravitational 
radius R g = GM/c 2 for masses in the 10 8 ~ 9 M Q 
is 1.5 x 10 13 ~ 14 cm. Taking a jet basis at 6 grav- 
itational radii one has typical dimensions in the 
range 2 x 10 14 ~ 15 cm. These results confirm that 
the blazar duty cycle is able of extremely rapid 
changes in flux both at 7 rays and optical frequen- 
cies and can have strong impact in the evaluation 
of the cosmic high-energy background (Giommi et 
al. 2006, Draper & Ballantyne 2009). 

Variability on time scale of minutes in the opti- 
cal band appears a rather accessible probe of the 
relativistic jet structure. In future observations 
of PKS 0537-441 and similar sources one should 
take into account that rapid variability of the type 
reported here appears a rare event and therefore 
very long monitoring should be scheduled. An im- 
portant improvement with respect to our photom- 
etry would be the actual simultaneity in all bands, 
which can be obtained by beam splitting, since one 
could probe the color of the short events. 
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Table 1: Example of photometry of PKS0537-441 from December 2004 to March 2009. The full version of 
the table is available online. 



MJD 


Exposure 


Filter 


Magnitude 


Error 


d 


s 




mag 


mag 


53357.17192 


600 


V 


16.20 


0.08 


55326.98610 


300 


V 


14.02 


0.11 


53360.16350 


30 


R 


14.93 


0.15 


55327.01266 


300 


R 


13.88 


0.03 


53358.18466 


600 


I 


15.28 


0.02 


55327.01645 


300 


I 


13.43 


0.06 


53676.17050 


30 


J 


13.04 


0.05 


55326.99834 


30 


J 


12.37 


0.03 


53677.29373 


30 


H 


12.62 


0.26 


55327.01354 


30 


H 


11.61 


0.05 


53453.02882 


30 


K 


10.71 


0.04 


55312.04763 


30 


K 


10.59 


0.05 



Table 2: Maximum and minimum intensity state of PKS0537-441. 



MJD 


filter 


Magnitude 


Flux(mJy) 


AMagnitudc 


55255.22365 


V 


12.95±0.09 


24.1 ±2.0 




53677.29965 


V 


16.87±0.25 


0.65 ±0.15 






V 






3.919±0.266 


55292.11168 


R 


12.08±0.07 


42.1 ±2.6 




53726.23140 


R 


16.80±0.22 


0.55 ±0.11 






R 






4.716 ±0.233 


53406.11486 


I 


12.00±0.08 


35.8 ±2.5 




53726.35053 


I 


16.37±0.24 


0.64 ±0.14 






I 






4.37 ±0.254 


55271.07350 


J 


10.78±0.12 


77.5 ±8.9 




53834.98987 


J 


17.78±0.25 


0.12 ±0.03 






J 






6.998 ±0.274 


54063.13561 


H 


10.86±0.03 


46.3 ± 1.4 




53835.04113 


H 


15.99±0.05 


0.41 ±0.02 






H 






5.129 ±0.062 


54063.14019 


K 


10.04±0.03 


64.5 ± 1.6 




53835.01757 


K 


15.37±0.06 


0.47 ±0.03 






K 






5.334 ±0.068 
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Tabic 3: A short event. 



MJD 


Filter 


Flux 


AT 


Tint 


Significance 






mjy 


min 


min 


a 


54038.29736 


J 


7.22±0.09 








54038.29944 


J 


7.69±0.06 










J 




3 


25.2±0.3 


5.4 
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